Automation of One-Loop Calculations with 
Go Sam : Present Status and Future Outlook* 

G. Cullen^ N. Greiner^'^, G. Heinrich^, G. Luisoni^, 
P. Mastrolia^, G. Ossola^'^ T. Reiter^ F. Tramontano^ 



^ School of Physics and Astronomy, 
The University of Edinburgh, Edinburgh EH9 3JZ, UK 

^ Max-Planck Insitut fiir Physik, 
Fohringer Ring, 6, D-80805 Miinchen, Germany 

^ Department of Physics, University of Illinois at Urbana-Champaign, 
1110 W Green Street, Urbana IL 61801, USA 

^ IPPP, University of Durham, Durham DHl 3LE, UK 

^ Physics Department, New York City College of Technology, 
City University of New York, 300 Jay Street, Brooklyn NY 11201, USA 

® The Graduate School and University Center, The City University of New York 
365 Fifth Avenue, New York NY 10016, USA 

^ Theory Group, Physics Department, CERN, 1211 Geneva 23, Switzerland 

Abstract 

In this presentation, we describe the GoSam (Golem/Samurai) 
framework for the automated computation of multi-particle scatter- 
ing amplitudes at the one-loop level. The amplitudes are generated 
analytically in terms of Feynman diagrams, and can be evaluated us- 
ing either D-dimensional integrand reduction or tensor decomposition. 
GoSam can be used to compute one-loop corrections to Standard 
Model (QCD and EW) processes, and it is ready to link generic model 
files for theories Beyond SM. We show the main features of GoSam 
through its application to several examples of different complexity. 

* Presented by G. Ossola at the XXXV International Conference of Theoretical Physics 
"Matter to the Deepest" : Recent Developments in Physics of Fundamental Interactions, 
Ustron 2011 
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1 Introduction 



The discovery potential of the experimental programs at the LHC relies 
heavily on the availability of higher order corrections for many relevant pro- 
cesses [l]. The searches for the Higgs boson and the compilation of related 
exclusion limits need precise calculations for Higgs' signal and background 
processes. Further, it will be very important to have precise theory pre- 
dictions at hand in order to constrain model parameters in the event that 
a signal of New Physics will be detected. Therefore, it is of major impor- 
tance to provide tools for next-to-leading order (NLO) predictions which are 
largely automated such that signal and background rates for a multitude of 
processes can be estimated reliably. 

Already some time ago, the idea of automating NLO calculations has been 
pursued with public programs like FeynArts [2] and QGRAF [3] for diagram 
generation and FormCalc/Loop Tools |4] and GRACE [5] for the automated 
calculation of NLO corrections, primarily in the electroweak sector. In spite 
of this important progress, until the last few years we did not observe a 
large production of calculations of one-loop amplitudes involving more than 
four external legs. Only very recently, conceptual and technical advances in 
multi-leg one-loop calculations allowed the calculation of six-point (6 22 and 



even seven-point 23,24 processes, and opened the door to the possibility 
of an automated generation and evaluation of multi-leg one-loop amplitudes, 
rather than creating a collection of hard-coded individual processes. 



Even if excellent process-specific programs are available, like MCFM 25 



27 and VBFNLO 28 , nevertheless it is desirable to have flexible tools at 



hand such that, in the same fashion already available at the tree-level 29 31 
any process which may turn out to be important can be promptly evaluated 
at NLO accuracy. 

Recently, we observed major advances in the direction of constructing 



packages for fully automated one-loop calculations, see e.g. 32 41 . Review- 
ing all the concepts that lead to these advances is beyond the scope of this 
presentatiorQ In the development of our computational tools, the OPP re- 
duction technique (43p4] and generalized D-dimensional unitarity fis] turned 
out to be the most crucial ingredients. 



The purpose of this talk is to present the program package GoSam 46 



which allows the automated calculation of one-loop amplitudes for multi- 



^ Additional information can be found in other talks presented in this conference 42 . 
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particle processes. The integrand is generated via Feynman diagrams, using 
QGRAF |3|, FORM |47 , spinney ES] and haggles H9\. The individual program 
tasks are managed by python scripts. The only task required from the user 
is the preparation of an "input card" in order to launch the generation of 
the source code and its compilation, without having to worry about internal 
details of the code generation. 

Concerning the reduction, the program offers the possibility to use ei- 
ther the Z)-dimensional extension of the OPP method, as implemented in 
SAMURAI 36 , or tensor reduction as implemented in Golem95C 50,51 



interfaced through tensorial reconstruction at the integrand level 52j, or a 
combination of both. 

GoSam can be used to generate and evaluate one-loop corrections in 
both QCD and electro-weak theory. Beyond the Standard Model theories 
can be interfaced using FeynRules 53 , 54 or LanHEP ^55j . The Binoth-Les 
Houches-interface [56j to programs providing the real radiation contributions 
is also included. 

In the following, we will provide a brief description of the main features 
of the code, with particular attention to the generation of the code and the 
various options to efficiently and automatically compute all rational terms. 
We will conclude the presentation with some examples of applications. 



2 Main features of GoSam 

GoSam produces in a fully automated way all the code required to per- 
form the calculation of one-loop amplitudes, by processing the information 
contained in an "input card" prepared by the user. The main steps in this 
process are: the generation of contributing diagrams, the optimization and 
algebraic manipulation to simplify their expressions, and the writing of a 
FORTRAN code ready to be used within a phase-space integration. The re- 
duction of unintegrated amplitudes to linear combinations of scalar (master) 
integrals is fully embedded in the process. 

In this section, we give a brief overview of some general operations per- 
formed by GoSam. A complete description of the framework, together with 
a detailed explanation of all features available in GoSam, can be found in 
Ref. Eel. 
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2.1 Diagram Generation 

For the diagram generation both at tree level and one-loop level we employ 
the program QGRAF js] . This program already offers several ways of excluding 
unwanted diagrams for example by requesting a certain number of propaga- 
tors or vertices of a certain type or by specifying topological properties such 
as the presence of tadpoles or on-shell propagators. Although QGRAF is a very 
reliable and fast generator we added another filter over diagrams by means 
of Python. This gives several advantages: first of all, the possibilities offered 
by QGRAF are not always sufficient to distinguish certain classes of diagrams; 
secondly, QGRAF cannot handle the sign for diagrams with Majorana fermions 
in a reliable way; and finally, in order to fully optimize the reduction, we want 
to classify and group diagrams according to the sets of their propagators. 
In our framework, QGRAF generates three sets of output files: an expres- 



sion for each diagram for FORM 47 , Python code for drawing each diagram 
and Python code for computing the properties of the diagram. The model 
information for QGRAF is either read from the built-in Standard Model file 



or is generated from a user defined LanHEP 55 or Universal FeynRules 
Output (UFO) [53] file. The Python program automatically performs several 
operations: diagrams whose color factor turns out to be zero are dropped; 
the number of propagators containing the loop momentum, the tensor rank 
and the kinematic invariants of the associated loop integral are computed; 
diagrams with a vanishing loop integral associated are detected and flagged 
for the diagram selection; all propagators and vertices are classified for the 
diagram selection; diagrams containing massive quark self-energy insertions 
or closed massless quark loops are specially flagged. 

Partitioning diagrams with similar structures and tracking their rank are 
very important operations in order to reduce the number of operations per- 
formed by the reduction and allow allow for a big gain in efficiency: after 
carrying out the tensor reduction for one diagram, all other diagrams that 
contain only a subset of the denominators are reduced with virtually no ad- 
ditional computational cost. This is true both in the OPP method f43| as 



implemented in Cut Tools 32 and SAMURAI [36j and in classical tensor 



reduction methods as implemented in Golem95C 50[ 51 , P JFRY 57\ 58 and 
LoopTools 



During this phase. Go Sam also generates a ETgX file with the drawings 
of all contributing diagrams. To achieve this task, we use our own implemen- 
tation of the algorithms described in Ref. (6O] and Axodraw IgT] to actually 
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draw the diagrams. 



2.2 Lorentz Algebra 

Concerning the algebraic operations performed by GoSam to render the in- 
tegral suitable for efficient numerical evaluation, one of the primary goals is 
to split the (4 — 2e) dimensional algebra into strictly four-dimensional ob- 
jects and symbols representing the higher-dimensional remainder. In GoSam 
we have implemented the 't Hooft-Veltman scheme (tHV) and dimensional 
reduction (DRED). In both schemes all external vectors (momenta and po- 
larisation vectors) are kept in four dimensions. Internal vectors, however, 
are kept in the n-dimensional vector space. We adopt the conventions used 
in 48 , where k denotes the four dimensional projection of an in general 
n dimensional vector k. The (n — 4) dimensional orthogonal projection is 
denoted as k. For the integration momentum q we introduce in addition the 
symbol /x^ = — g^, such that 

q=q+q=q-^. (1) 
We also introduce suitable projectors by splitting the metric tensor 

g,u^^,u^~,u^ ^^.^^^^0, g>; = A, r,=n-A. (2) 

GoSam contains a library of representations of wave functions and prop- 
agators up to spin two. The exact form of the interaction vertices is taken 
from the model files. 

Once all wave functions and propagators have been substituted by the 
above definitions and all vertices have been replaced by their corresponding 
expressions from the model file, all vector-like quantities and all metric ten- 
sors are split into their four-dimensional and their orthogonal part. As we 
use the 't Hooft algebra, 75 is defined as a purely four-dimensional object, 
75 = i^^iupal^l^l^l" ■ By applying the usual anti-commutation relation for 
Dirac matrices we can separate the four-dimensional and (n — 4)-dimensional 
parts of Dirac traces. 

While the (n — 4)-dimensional traces are reduced completely to products 
of (ra — 4)-dimensional metric tensors g^'^ , the four-dimensional part is treated 
such that the number of terms in the resulting expression is kept as small as 
possible. Any spinor line or trace is broken up at any position where a light- 
like vector appears. Furthermore, Chisholm identities are used to resolve 
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Lorentz contractions between both Dirac traces and open spinor lines. If any 
traces remain we use the buih-in trace algorithm of FORM 47 . 



2.3 Treatment of R2 terms 



In the numerator of a one-loop diagram, terms containing the symbols fi or 



e can lead to a so-called R2 term 62 . Therefore the numerator function can 
be written as, 

Af{q,fi^,e) =A/'o(g,Ai^)+eM(g,/i2)+£W2(g,/i^). (3) 

It is useful to observe that the terms A/i and A/2 in Eq. ([s]) do not arise 
in DRED, where only terms containing /i^ contribute to i?2- Instead of 



we 



relying on the construction of R2 from specialized Feynman rules 63 66 
can generate the R2 part along with all other contribution using automated 
algebraic manipulations. 

The code offers the option between the implicit and explicit construction 
of the i?2 terms. The implicit construction uses the splitting of Eq. ^ and 
treats all numerator functions A/i on equal grounds. Each term in Eq. ^ 
is reduced separately and the results are added up taking into account the 
powers of e. The explicit construction of R2 is based on the fact that the non 
purely 4-dimensional part of the numerator function contains powers of /i^ 
or e, and the expressions for the corresponding integrals are relatively sim- 
ple and known explicitly. Therefore, after separating it using the algebraic 
manipulation described before, the (n — 4) dimensional part is computed ana- 
lytically whereas the purely four-dimensional part is passed to the numerical 
reduction. This approach also allows for an efficient calculation of the R2 
alone. 



2.4 Reduction to scalar (master) integrals 

GoSam allows to choose at run-time (namely without regenerating the code) 
the preferred method of reduction. Available options include the integral- 
level D- dimensional reduction, as implemented in SAMURAI, or traditional 
tensor reduction as implemented in Golem95C interfaced through tensorial 
reconstruction at the integrand level, or a combination of both. 

Concerning the scalar (tensorial) integrals [67, 68 , GoSam allows to 
choose among a variety of options, including QCDLoop 69 70 , OneLoop 71 
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Golem95C |50||5T), plus the recently added PJFRY |57|(58| and LoopTools |4||59 
Among these codes, OneLoop and Goleni95C also fully support complex 
masses. 

For details about the reduction methods, we refer the reader to previous 



presentations 37 39 or the original articles. 



3 Examples 

The GoSam codes have been tested on several processes, starting with QCD 
2 — )■ 2 NLO amplitudes, but also on more challenging 2 — )■ 4 (not counting 
decays) in the final state. Some examples are depicted in Table 1. The full 
list of processes, with the details of all comparisons performed, is given in 
Ref. 



46 . 



Process 


Checked with Ref. 


ud — i- e^Ue g 


40 






e+e- e+e-7 (QED) 


72 






pp ^ Hti 


40 






pp W^W^jj 


20 






pp bbbb 


15 


16 


34 




pp W^W^bb 


34 


40 






uu —J- tibb 


34 


40 




gg — tibb 


34 


40 




ud W^ggg 


34 







Table 1: Some of the processes computed and checked with GoSam 



3.1 BLHA interface, GoSam, and SHERPA 

The BLHA interface allows to link GoSam to a general Monte Carlo event 
generator, which is responsible for supplying the missing ingredients for a 
complete NLO calculation of a physical cross section. Among those, SHERPA 73 
offers the possibility to compute the LO cross section and the real corrections 
with both the subtraction terms and the corresponding integrated counter- 



parts |74] - [76] . Using the BLHA interface, we linked GoSam with SHERPA 
to compute physical cross section for + 1-jet at NLO. 
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We tested our results producing distributions for inclusive and exclusive 
p_L and f] of the jet, Ht, and for p± and rj of the leptons (details can be found 
in Ref. [46]). All distributions are in agreement with the ones produced using 
SHERPA in combination with MCFM. 



1 1 1 1 , 1 1 1 1 , 1 1 1 1 , 1 1 1 


1,1111,111 


1 1 1 1 _ 
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Figure 1: Comparison of the inclusive transverse momentum of W~ + 1 jet 
between Sherpa interfaced with GoSam and Sherpa interfaced with MCFM. 
For comparison purpose we chose the /cT^-algorithm with pt^min = 50 GeV. 
The bars indicate the statistical Monte Carlo error. 



3.2 GoSam and Neutralino Pair Production 

As an example of the usage of GoSam with a model file different from the 
Standard Model, we calculated the QCD corrections to neutralino pair pro- 
duction in the MSSM. A calculation of the total cross sections for neutralino 



pair production at the LHC is also presented in Ref. 77 . The model file has 
been imported via the UFO interface. To import such files within the GoSam 
setup, all the user has to do is to give the path to the corresponding model 
file in the input card. 

In this example, we combined the one-loop amplitude with the real ra- 
diation corrections to obtain results for differential cross sections. For the 
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infrared subtraction terms we employed MadDipole [78|[79], while the real 
emission part is calculated using MadGraph/MadEvent |80] . The virtual 
matrix element is renormalized in the MS scheme, while massive particles 
are treated in the on-shell scheme. The renormalization terms specific to the 
massive MSSM particles have been added manually. For the SUSY parame- 
ters we use the modified benchmarks point SPSlamod suggested in 81 , and 
use a/s = 7TeV. 



^0^0 invariant 

A.1 A.1 



In Fig|2] we show the differential cross section for the m 
mass, where we employed a jet veto to suppress large contributions from the 
channel qg — )■ XiXiQ which opens up at order a^a^, but for large p^^* belongs 
to the distinct process of neutralino pair plus one hard jet production at 
leading order. 




200 400 600 



1200 1400 1600 
m(N11,N12) [GeV] 



Figure 2: Comparison of the NLO and LO fTixIxl distributions for the process 
PP ~^ XiXi with a jet veto on jets with p^^* > 20 GeV and rj < 4.5. The band 
gives the dependence of the result on fi = fip = fi^ between /io/2 and 2/xo. 
We choose fiQ = irtz- The black line gives the bin error for the value at the 
central scale. 
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4 Outlook and Conclusions 



In the last five years, we observed major advances in our understanding of 
one-loop scattering amplitudes. Aside from improvement on standard tenso- 
rial techniques, the development of unitarity-based approaches, paired with 
the decomposition at the integrand level contained in OPP method, changed 
the landscape of this field, favoring the calculation of NLO amplitudes for 
several challenging processes and the development of new theoretical frame- 
works and tools for such calculations. 

For quite a long time, tree-level calculation have been fully automated 



and included in flexible multi-process tools 29 , 30 . The level of automation 
achieved by one-loop calculations is suggesting the possibility of a similar 
success also for the NLO. One of the natural hopes for the future is to de- 
vise ways of extending and generalizing what we understood so far about 
one-loop amplitudes, in order to develop tools and methods for higher-order 



calculations 82 -84 , but plenty of work is still needed to achieve this goal. 

In this presentation, we illustrated the main features of GoSam, a new 
program package for the fully automated evaluation of one-loop scattering 
amplitudes in any renormalizable quantum field theory. In its present form, 
GoSam can be used to calculate one-loop corrections both in QCD and 
electro-weak theory and offers the flexibility to link general model files for 
theories Beyond the Standard Model. The amplitudes are generated in terms 
of Feynman diagrams and the reduction to master (scalar) integrals can be 
performed in several ways, which can be selected at run-time. 

We presented several examples of one-loop calculations performed within 
the GoSam framework, as well as preliminary results of interesting appli- 
cations such as the interface with SHERPA. These examples demonstrate 
the great flexibility, together with a competitive timing, of GoSam. We are 
looking forward to tackle more challenging calculations and interfacing with 
other existing tools in the coming months. 
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